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Abstract ; The  2 fol d mul t i pi ic i ty  of  solutions  found  by  Bo  ley 
& Golub  in  their  study  of  the  inverse  problem  for  N x N symmetric, 
pentadiagonal  matrices  contrasts  with  the  unicity  of  the  solution 
of  the  inverse  problem  for  an  inhomogeneous,  discrete  beam.  The 
reason  for  this  discrepancy  is  elucidated  and  can  be  traced  to  the 
different  properties  of  the  spectral  data  used  in  the  two  cases. 
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1.  Introduction.  Recently,  Boley  S Golub  [1]  have  examined  the  problem  of  re- 
constructing an  NxN,  symmetric,  pentad iagonal  matrix  from  its  spectral  data. 

The  formulation  of  this  inverse  problem  follows  the  line  pioneered  by  Hochstadt 
[2]  for  tridiagonal  matrices.  Namely,  the  given  spectral  data  are  made  up  by 
the  eigenvalues  of  the  original  pentad iagonal  matrix,  as  well  as  the  eigenvalues 
of  two  closely  related  matrices  obtained  by  deleting  respectively  the  first  one 
and  two  rows  and  columns.  Whereas  the  equivalent  procedure  yielded  a unique 
solution  for  tridiagonal  matrices  [2],  Boley  6 Golub  found  a 2^  ' multiplicity 
of  solutions  for  the  pentad iagona 1 case. 

This  result  might  seem  puzzling  at  first.  Indeed,  pentadlagonal  matrices 
can  be  viewed  as  finite  difference  analogues  of  certain  iith  order  differential 
operators,  such  as  that  governing  the  vibration  of  a beam.  And,  in  my  investi- 
gation of  the  Inverse  problem  for  the  discrete  beam  1.3],  I showed  that  given 
three  so-called  sympathetic  spectra,  the  solution  of  the  inverse  problem,  if 
it  existed,  was  unique. 

The  aim  of  the  present  paper  is  to  clarify  this  apparent  contradiction 
between  Boley  S Golub's  result  and  mine.  The  key  to  the  paradox  lies  in  the 
fact  that  Boley  6 Golub  used  three  spectra  which  are  not  sympathetic,  i.e.  which 
are  not  equivalent  in  information  to  the  data  contained  in  the  impulse  response. 

In  order  to  illustrate  this  important  point  about  the  need  for  the  spectral 
trio  to  be  interrelated,  1 shall  consider  a specific  inverse  problem  for  the 
vibrating  beam  with  three  non-sympathetic  spectra.  More  specifically,  the 
spectral  data  will  consist  of  the  natural  frequencies  of  vibration  of  the  beam 
in  the  following  three  configurations:  (i)  clamped-c lamped,  (ii)  clamped-supported 
and  (iii)  clamped-free.  Eminently  reasonable  though  this  choice  of  spectra  may  be 
from  the  engineering  point  of  view,  it  has  the  drawback  that  these  three  spectra 
are  not  sympathetic  and  cannot  insure  a unique  solution.  In  fact,  for  the  Nth 
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discretized  version  of  this  problem,  we  shall  see  that  there  are  2^  ^ solutions, 
thus  arriving  at  a multiplicity  of  solutions  similar  to  that  found  by  Boley  6 
Golub. 

The  outline  of  the  paper  is  as  follows.  In  §2,  we  shall  derive  the  basic 
equations  for  a d i screte  beam:  this  simple  mechanical  system  will  provide  valu- 
able Insights  into  the  inverse  problem.  §3  contains  all  the  necessary  ingredients 
for  a consideration  of  the  inverse  problem.  Finally,  in  %k  we  consider  the  partic- 
ular inverse  problem  for  which  the  spectral  data  are  associated  with  the  natural 
frequencies  of  vibrations  in  the  clamped/clamped,  clamped/supported  and  clamped/ 

free  configurations.  We  shall  show  that  with  such  data,  the  inverse  problem  has  a 
N-i 

2 fold  multiplicity  of  solutions. 
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2.  The  dtscrete  beam.  Consider  a beam  of  length  L characterized  by  a variable 
flexural  rigidity  El  and  a variable  density  p.  If  y(x,t)  stands  for  the  dis- 
placement of  the  central  line,  then  the  classical  equation  governing  the  infini- 
tesimal oscillations  of  this  beam  is 

" El— , 0 < X < L.  (1) 

3x^  3x^  3t 

For  the  sake  of  definiteness,  we  shall  assume  that  the  left  end  of  the  beam  is 
clamped,  i.e. 


y(o.t)  - y (0,t) 


(2) 


Let  us  now  consider  the  following  conceptual  experiment:  at  time  t ■ 0, 
we  strike  impulsively  the  right  end  of  the  beam  which  is  assumed  to  be  stress- 
free.  We  then  monitor  the  ensuing  displacement  as  well  as  the  slope  of  the 
central  line  of  this  right  end.  To  carry  out  this  program  mathematically,  we 
must  solve  equation  (1)  subject  to  (2)  and  to 


y(x,0)  - ^ y(x,0)  * 0, 

2 2 

~ y(L,t)  * ^ El  y(L,t)  - F5(t)  - 0, 


(3) 


3x‘ 


3x 


3x 


in  order  to  deduce  the  impulse  response,  viz.  y(L,t)  and^^(L,t).  This  is  best 
accomplished  by  working  In  the  frequency  domain;  to  that  effect,  we  define 


u (x,a)) 


e'‘^^  y(x,t)  dt. 
o 


More  importantly,  in  preparation  for  the  discretization,  it  is  convenient 
to  transform  (i)  into  a system  of  first  order  equations.  Therefore,  let  us  define 
0(x,a))  to  be  the  slope  of  central  line,  i.e. 

u'  - 9 (5) 

and  T (x,aj)  to  be  the  moment  about  the  central  line,  i.e. 

9'  ■ f{x)T. 


(6) 
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In  the  above  formulas,  a prime  denotes  differentiation  with  respect  to  x and 
f(x)  ■ 1/El  is  a measure  of  the  "limpness"  at  x.  To  complete  our  transformation 
to  a first  order  system  of  equations,  we  introduce  (|)(x,oj)  which  represents  the 


force  applied  to  the  central  line.  Consequently, 

t'  = <J)  (7) 

With  these  new  variables,  Newton's  law  reads 

(j)'  = -(i)^pu.  (8) 

The  boundary  conditions  (2)-(3)  can  now  be  written  thus: 

u(0)  “ 9(0)  » 0,  (9) 

t(L)  - (J)(L)-F  = 0.  (10) 


We  shall  be  concerned  with  the  discrete  version  of  the  system  (5)“(8) 
subject  to  conditions  (9)“(10).  The  simplest  and  most  physically  meaningful 
way  to  discretize  this  problem  consists  in  restricting  p(x)  and  f(x)  to  be  sums 


of  delta  functions,  viz. 

N 

p(x)  • I m.  6(x-x.) 
i-1 

N 

f(x)  - f.  5(x-x.) 
i-1  ' ' 


(11) 


segments  are  connected  by  joints  of  mass  mj  and  concentrated  limpness  fj  located 
at  X - X..  For  instance,  we  can  think  of  these  joins  as  clothespin-like  devices. 
We  can  see  from  (6)  and  (8)  that  in  the  interval  (xj_j»  ■’^j)*  l^^li  ® ^nd  i|) 


are  constants,  i .e. 
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e(x,u))  - 0._j 

(p(x,lj))  - 


for  xe (x._^ , X. ) 


Substituting  these  expressions  in  (5)  and  (7)  we  deduce  that 
u(x,a))  * u.^^  + (x-x._j)  0._j 
t(x,U>)  - T._j  - (X-X._j) 

We  get  the  discrete  equations  by  examining  the  jump  conditions  at  x.< 
seen  from  (5)”(8),  these  conditions  are: 

[u]j  - 0 

[0].  * f.  T(x.,a)) 

[r],  - 0 

2 

[(})]  J * -w  mj  u{xj  ,0)) 

where  [u].  stands  for  u(x.+0,(i))  - u(x.-0,a)),  etc. 


Replacing  (14)-(15) 


nto  (16)  yields 


u . 

1 

u. 

1 - 

1 ^ 

*'i-l®i-l 

0. 

1 

- 

0. 

1 - 

1 ^ 

f.T. 

1 1 

T. 

1 

- 

T._ 

1 

^i-l'^i-1 

<P, 

- 

4>i. 

1 

a)^m . u . 

1 1 

Since  u and  T are  continuous,  it  is  clear  that 
Uj  - u(x.,u))  , 

Tj  - T(xj,a)), 

whereas  by  the  very  definition  in  (lA) 

0.  - 0(x.+O,a)), 

4».  - <))(xj+0,cj) . 

To  complete  the  formulation  of  the  discrete  d i rect  problem,  we  shall 

following  boundary  conditions: 

u *0  *0, 

o o ’ 


‘n  ■ Vf 


(1M 

(15) 

As  can  be 

(16) 

(17) 


impose  the 

(18) 


(19) 
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Before  proceeding  with  the  computation  of  the  impulse  response,  viz.  Uj^  and  0|^, 
it  is  perhaps  worth  writing  (17)  as  a single  forth  order  difference  equation, 
name  1 y 

1 ..  r 1 ] . I—  ! ; .. 


1 1 1 1 1 ^ 

! o“T  ^i+2~^f  Si  i * f Si  ^ ^ f 9 2 ^ f £ £ i + 1 

i+ri*'i+i  ' ’^i+ri  i+1  ’^i+ri  ’^i^i  ^i^i-ri 


I J 2 1 1 

0 o~Z  * f 0 Z * f Si  £ ''’f  £ 2 f Si  2^ 

i+ri  ^i  i ’^i^i-ri  i*i-i  ^i-ri-i 


1 1 1 1 ' 
"■ffo  0 * f Si  ^ f £ f £ £ "^  *^i”l 


f £ £ “{-2  “ '"i‘'i* 

^-ri-2^i-l  ' ^ ' ' 


Now,  if  we  introduce  the  field 


V,  =>  m.  u. 
t I I 


then  (18)  can  be  written  as 


C..2  Vi-z  + ^-l''i-l  ^''i  ^''i+l  ‘^i' 


whe  re 

1 

r > 

+ ' + 

2 ^ 

+ 

1 

o • 

1 m. 

1 

f £ 2 
, i + 1 1 

f.£i2 

^ 1 

f.£.  2 

1 1-1 

f.  ,£. 

1-1  1-1 

f 

-1  1 

1 

V 

+ 1 

+ ’■  1 

' /-  5 

(mjmj+i) 

f £ £ 
Ji+l^l*i+l 

f..,£.2 
1+1  1 

f.^,2 

1 1 

^i^i-l^i 

9 

{21a) 


(21b) 


/m.m...\^  f •.|X».£,| 

( I 1+2)  1 + 1 I 1 + 1 


(21c) 


Thus,  the  inverse  problem  for  a discrete  beam  corresponds  to  an  inverse  problem 
for  a symmetric  pentad lagona 1 matrix  in  which  the  elements  of  the  first  off- 
diagonals  are  negative,  all  other  elements  being  positive. 
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One  final  remark  regarding  the  discretization.  No  particular  significance 
should  be  attached  to  the  choice  of  the  clamped  boundary  conditions  (18).  How- 
ever, having  made  that  choice,  it  is  easy  to  see  from  (17)  that 

u,  - 0,  (22) 

i.e.  the  first  joint  does  not  move!  As  a result,  the  mass  m^  does  not  enter 
into  the  problem.  In  fact,  for  this  clamped/free  case,  the  eigenvalue  problem 
for  the  pentad iagonal  matrix  looks  as  follows: 


(23) 


where 


'n-1 


"N-l 


Vl*'N-l^ 


9 9 

•rN-2  N-l 


2,  2 

•1  N-2 


f 2 2 

N-2  N-2 


-1 


"n-1 


^N-l^N-1^ 


^N-1^N-2^N-1 


"’n^n-i*'n-i^ 


We  shall  not  pursue  the  problem  via  the  pentadiagonal  matrix  since  this  tends 
to  shroud  the  physics. 

3.  The  wedge  product.  The  direct  probiem  consists  in  finding  Uj^  and  9|^  given 

N N N 

the  beam  characteristics  solution  to  this  problem  is 

straightforward.  We  introduce  two  linearly  independent  fundamental  solutions 


I 
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of  the  difference  equations  (17).  Using  a superscript  to  label  these  solutions, 


we  define  [uf’^  0.^’^  (()p^]  and  [up\  9p\  Tp\  or 

succintly  u^^  and  u^^^  in  such  a way  that 

- [0  0 I 0]^  , 


up^  = [0  0 0 1]''’, 


where  T stands  for  transpose. 

Then  . 

u.  = au:  ^ + bu:  ^ 

~i  ~i  ~i 


more 


(2i.) 


(25) 


In  order  to  evaluate  the  constants  a and  b,  we  make  use  of  the  boundary  conditions 

(19)  at  X,,,  namely 
N 


(26) 


,(2)- 

N 

a 

3 

0 

‘*’n 

*(2) 

‘'’n 

_ 

b 

F 

Omitting  the  trivial  intermediary  calculations  we  find  that  the  constituents  of 
the  impulse  response  can  be  expressed  thus: 


and 


0. 


-F 


>F 


{^N’  '^N} 


(27a) 


(27b) 


where  from  now  on  we  shall  use  a bracket  as  a shorthand  notation  tor  determinants 


of  certain  2x2  matrices,  namely 

{,.  b}  - ,<’>  b'^>  - ,(2)  ,0) 


(28) 
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Bracket  expressions  will  be  formed  for  each  point  x. , e.g. 


r 1 (1)  (2)  (2)  (1) 

{u,,  T.}  » u'  ' T.  ' - u.  ' T.; 


I ’ I 


These  bracket  expressions  are  related  to  the  so-called  wedge  product  of  Uj 

(2) 

and  u.  . Given  the  4x2  matrix 

~ I 


(1) 


u.  = [u!’^  , 

~l  ~l  * 


the  wedge  product 


Z.  . u<')  A 

~ I ~ I - I 


(29) 


is  a vector  with  components  . (a  = 1,  2,  . . .,6)  equal  to  the  determinants 
of  the  matrices  obtained  from  U.  by  deleting  two  rows.  If  the  label  a is 
attached  to  the  chosen  rows  in  a lexicographic  manner,  then  it  is  clear  that 


•l,i 


{u 


^2,1  = 


{u 

{0 


'3,  i 

■^.i 


Z.  . = {9 

5 f * 

^6.i 


. 9,} 

. T.} 

f <J>j} 
. T,} 

» 'I’j  } 


(30) 


Rather  than  working  with  u.^^^  and  u!^^  whose  components  are  polynomials 

2 

in  (i)  and  then  forming  the  bracket  expressions,  a risky  numerical  procedure 
in  view  of  the  cancellations,  it  is  preferable  to  work  with  the  brackets  ab  initio. 
Our  next  task  then  is  to  derive  equations  for  Z,  .,...,  Z,  ..  Using 

1,1  b, I 

the  difference  equations  (17),  it  is  a simple  matter  to  deduce  that 

^1,1  " ^1,1-1  ^i^2,i  ’ 


Z,  , - Z-  . , - a,  ,Z,  . , + S,.  ,Z,  . , - 7 Z_  . , 

2,1  2,1-1  i-1  3*i“l  I"!  4,1-1  1-1  5.i“l 


(31b) 


10  - 


z 

z 

z 

z 


3.f 

^.5 

5.1 

6. 1 


= Z 

- Z 

= Z 

= Z 


3.1- 1 

4.1- 1 

5.1- 1 

6. 1- 1 


^1-1^5, 1-1 


0 7 

^i-r5,i-i 


+ f.Z,  . , + w m.Z,  . , 

I 6, I - 1 I 1 , 1 


+ oj  m.Z«  . • 

I 2,1 


(31c) 

(31d) 

(31e) 

(31f) 


From  the  boundary  conditions  (24),  It  follows  that 


Z = '5  , 
a,o  ao 


ot  = 1 , 2,  . . • , 6 


(32) 


The  above  equations  can  be  simplified  slightly. 
(31d),  we  see  that 

^3,1  ^4,1  “ ^3,1-1  ^4,1-1 


which,  In  view  of  (32),  Implies  that 

Z,  . = -Z,  . . 

3,1  4,1 


We  can  therefore  eliminate  Z.  . and  rewrite  equations 

3,  I 

Z,  • = Z,  . , + f.Z«  . , 

1,1  1,1-1  <2,1  ’ 


Z,  . = Z-  . , + 22..  , Z.  . , - J..  , Zr  ■ , 
2,1  2,1-1  1-1  4, 1-1  1-1  5,1-1 


■li  i “ j 1 - *>!  i^c  • I 

“t,!-!  1-1  5,1-1, 


Z,.  , * Z_  , , + f.Z,.  . , + 0)  m.Z,  . 


■5.1  5,1-1  ^ 'r6,l-l 


• 4*  « 

I 1 , 


Z/;  : * 2,  . , + 0)  m.Z-  . . 

6,1  6,1-1  <2,1 


Indeed,  from  (31c)  and 


(33) 

(31)  as 

(34a) 

, (34b) 

(34c) 

(34d) 

(34e) 


These  are  the  equations  which  are  best  suited  for  solving  the  inverse  problem. 
It  is  important  to  realize  that  the  components  of  Z^  . are  not  Independent. 
Indeed,  in  addition  to  the  linear  identity  (33),  these  components  satisfy  the 
following  quadratic  identity: 


(35) 


z,  .z,  . - z.  : - z-  ,Zc  . = 0 

1 , 1 6,  I 4,1  2, 1 5, 1 


This  identity  can  be  derived  from  the  equations  (3^).  However,  such  a derivation 

hides  the  fact  that  this  identity  is  algebraic  in  nature.  It  has  a long  history 

which  goes  back  to  work  done  in  the  eighteenth  century  on  "vanishing  aggregates 

of  determinants"  (see  e.g.  [4],  bottom  of  p.  50).  This  quadratic  identity  will 

play  a crucial  role  in  whether  or  not  the  inverse  problem  has  a unique  solution. 

We  close  this  paragraph  by  returning  to  the  determination  of  the  impulse 

response.  In  view  of  (27),  this  is  tantamount  to  the  determination  of 
2 2 2 

^2  ^4  ^ ^6  N Starting  with  the  conditions  (32),  we  can 

2 2 

find  these  quantities,  as  well  as  Z^  (oj  ) and  Z^  (u)  ),  by  solving  the  difference 

2 

equations  (34).  For  N>2,  these  quantities  are  polynomials  in  of  the  following 


Z,  ) = Z,  „(0) 


(1  - ^) 
n 


(36a) 


Z.,  „(w  ) - Z, 


(1  - ^ ) , 

X ^ 
n 


(36b) 


TT  (1  -isi)  . 

1 1 


(36c) 


(!-%■)  . 


(36d) 


(1  - ) . 


(36e) 


The  zeros  of  these  polynomials  are  the  squares  of  the  eigenfrequencies 
of  the  discrete  beam  in  various  vibrating  configurations.  These  vibrating  con- 


figurations are  such  that  the  left  end  is  always  clamped,  whereas  the  conditions 
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at  the  right  end  differ.  For  instance  {k  {X  }*l!  ^ and  (C  }*!!  ' correspond 

n I n I n I 

respectively  to  the  cases  In  which  the  right  end  is  clamped,  supported  and 

stress  free.  The  constants  ’ • *»  ^6  appearing  in  (36)  are 

2 

easily  deduced  by  setting  o)  = 0 in  (3^)  and  solving  the  resulting  difference 
equations.  This  calculation  shows  that 


Zi^^(O)  » Q=  - f.  I (x.'Xj)^  fj,  (N>2) 
i=2  j=I 


(37a) 


“ -'’2^-V  (V-;)'  . (N>2) 

1 = 1 


(37b) 


N-I 

* ■'’1= 


(37c) 


1 . 


(37e) 


4.  The  inverse  problem.  As  previously  mentioned,  in  view  of  our  particular 
choice  of  boundary  conditions  at  the  left  end,  m^  does  not  enter  into  the 
problem.  Therefore,  the  solution  of  the  inverse  problem  consists  in  recovering 
the  3N-1  beam  characteristics  ’ ^^i^O  ' * 

In  order  to  arrive  at  a unique  solution,  we  must  use  impulse  response 
data  together  with  two  additional  constraints  [3].  We  have  seen  that  the 

impulse  response  data  are  equivalent  to  the  three  sympathetic  spectra 

,N-2  t ,N-2  ,N-1  j , 

^^n  I ’ ^^n  ] ’ ^^n  ) 


N-I 

Pi  “ I (x,,-x.)  f.  , 

1 .‘'i  N I I 

1=1 

N-1  2 

Py  " I.  ^i 


(38) 


i-1 


(39) 
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To  these  3N-3  pieces  of  information,  we  add  another  sum,  say 
N-1 


•o-.l 


f. 


(1*0) 


i-1 

and  a normalization 
N-I 


«L  - I ^i“'-  • 
^ i*0 


(41) 


Of  course,  we  could  translate  these  requirements  into  analogous  ones  for 

pentadiagonal  matrices  of  the  type  (23).  But,  without  the  underlying  connection 

with  the  beam,  these  conditions  would  not  be  very  meaningful. 

We  shall  not  review  the  procedure  necessary  for  solving  the  inverse 

problem.  Suffice  it  to  say  that  this  procedure  consists  of  two  steps  [3]. 

2 

The  first  step  is  the  evaluation  of  the  missing  polynomials  Zj  (oj  ) and 
2 

Zj  Hj(<*'  ) ''■3  quadratic  identity  (35).  These  polynomials  can  be  determined 
uniquely  provided  that  the  given  data  are  associated  with  the  impulse  response. 
Indeed,  substituting 

N-2  , 

(42) 


7 ( 2\  2< 
<-0 


and 


• "o  * i!  V'" 


(43) 


K"! 


into  the  quadratic  identity  and  equating  the  coefficients  of  various  powers 
2 

of  0)  equal  to  zero  we  can  derive  a set  of  2N-2  1 inear  equations  for  the  un- 

N*  1 N*  2 2 2 

knowns  and  (q|^)Q  • Having  determined  Zj  nj(<^^)  and  Z.  |ij(w^)  and 

2 

thus  knowing  all  the  components  of  Z ^ (oj  ) we  can  proceed  with  the  second 

Ct  f N 

2 

step.  This  consists  in  finding  m^^,  fj^,  and  Z^  n-1  ^ from  (34)  by 

means  of  suitable  operations  with  polynomials  such  as  divisions  and  additions. 
Once  again,  the  results  of  these  operations  are  unique.  By  repeating  the 


- li,  _ 

same  operations  for  the  successive  intervals  the  unique  solution  to  inverse 
problem  is  obtained. 

The  above  mentioned  procedure  breaks  down  if  the  given  spectra  are  not 

sympathetic.  In  order  to  illustrate  the  changes  brought  about  by  data  which 

are  not  related  to  the  impulse  response,  let  us  consider  the  case  in  which 
2 2 2 

Z,  ,,(a)  ),  Z-  ,,(aj  ) and  Z,  ..(co  ) are  given.  Referring  to  (30)  for  the  defini- 

i|N  DyN 

tions  of  Z,  ,,,  Z»  ,,  and  Z,  it  is  clear  that  the  zeros  of  these  polynomials 

correspond  to  the  natural  frequencies  of  vibrations  of  the  beam  in  the  clamped/ 

clamped,  clamped/supported  and  clamped/free  configurations.  Following  the 

2 2 

same  procedure  as  before,  we  first  attempt  to  find  Z^  ) and  Z^  n ^ 

which  we  write  thus: 

2 2< 

Once  again,  the  quadratic  identity  provides  us  with  2N-2  equations  for  the  un- 

N“  2 N*  2 

knowns  at  hand,  namely  {q  and  (s  }_  . However,  these  equations  are 

no  longer  linear  equations  for  q^  and  s^  but  rather  linear  for  q^  and  quadratic 

for  s 1 By  eliminating  the  q 's  we  can  deduce  N-1  quadratic  equations  for  the 
1C  1C 

S|^.  Appealing  to  a theorem  of  Bezout^  (see  e.g.  [5],  p.  10)  there  are  in 

general  2*^  ' solutions  to  these  equations  i.e.  2*^  ' polynomials  Z^ 

2 N*  1 

Zj  l^(a)  ) compatible  with  the  given  data.  For  each  one  of  these  2 pairs,  we 

2 

can  find  corresponding  values  of  m^^,  f^^  and  Z^  n-1  ^ means  of  the 

N- 1 

usual  polynomial  manipulations.  Thus,  we  generate  2 solutions  to  the  inverse 
problem. 

The  second  definition  of  sympathetic  spectra  [3]  is  now  quite  natural. 

Such  spectra  are  the  zeros  of  those  three  polynomials  for  which  the  remaining 

^Bezout's  theorem:  N polynomial  equations  of  degrees  n^  n.,  . . • » 
in  N variables  have  in  general  n^  n^  . • • nj^  common  solutions,  when  the  number 

is  greater,  it  is  Infinite. 
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two  polynomials  can  be  uniquely  determined  by  the  quadratic  identity.  For 

2 2 2 

instance,  the  zeros  of  2j  ),  (to  ) and  ^2  ^ ^ sympa- 

thetic spectra.  In  their  investigation  of  the  inverse  problem  for  symmetric, 
pentadiagonai  matrices,  Boley  6 Golub  chose  their  spectral  data  primarily 
for  convenience  and  without  regard  to  whether  or  not  the  three  spectra  were 
sympathetic.  As  a result,  the  fact  that  they  found  a multiplicity  of  solutions 
is  no  longer  surprising. 
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